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We first consider the origin of colour deconfinement, comparing in particular sponta- 
neous symmetry breaking and percolation. Next we specify the onset of deconfinement 
^ ■ in nuclear collisions through parton percolation, and with parameters determined at SPS 
energy, we give quantitative predictions for charmonium suppression at RHIC. Finally we 
show that the parton cascade model leads to similar predictions, starting from a micro- 
OO I scopic description of the space-time evolution of the collision. 

5 : 1. INTRODUCTION 
^ , 

Q-i! The interesting new results on the threshold and the further pattern of anomalous 
{iLiI J/ip suppression |jl],0] make it appropriate to take stock of what we can say about the 
onset of deconfinement in nuclear collisions and what we can predict for the forthcoming 
RHIC experiments. The aim of this survey is therefore to formulate the conditions for 
. deconfinement, use the SPS results to determine the crucial onset parameters, and then 
I calculate what should be seen at RHIC. 

I Finite temperature QCD predicts that strongly interacting matter will undergo a tran- 
sition at T = Tc] at this point, hadronic matter, as the colour insulator phase of QCD, 
turns into the colour-conducting quark-gluon plasma. To throw some light on how and 
why this happens, we begin with a study of physical patterns for deconfinement. 

1.1. Spontaneous Symmetry Breaking vs. Cluster Percolation 

For an infinite bare quark mass, the QCD Lagrangian leads to pure colour SU{N) 
gauge theory. Here the expectation value of the Polyakov loop L{T) provides an order 
parameter suitable to define the confinement- deconfinement phase transition, driven by 
the underlying center Zn C SU{N) symmetry of the Lagrangian Confined states 

share this symmetry, but it is spontaneously broken in the deconfined phase. Spin systems, 
in particular the iV-state Potts model {N=2 is the Ising model), show the same kind of 
critical behaviour, based on the spontaneous breaking of a global Z^- symmetry, with the 
'magnetisation' m{T) as order parameter. This similarity has led to the conjecture that 
the two systems are in fact in the same universality class of critical behaviour which 
is indeed nicely confirmed in lattice studies B. 



*DifFerent parts of the material presented here are based on joint work with S. Fortunate, H.-W. Huang, 
D. Kharzeev, M. Nardi and D. Srivastava. 
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To arrive at an alternative physical definition of deconfinement, we will make use of the 
relation between spin and gauge systems. The critical behaviour of the Ising model can 
today be characterized either as 

• spontaneous magnetisation, disorder/order transition connected to the 
spontaneous Z2 symmetry breaking for T < Tc, or as 

• percolation of clusters of parallel spins joined by suitable bonds, i.e., in terms of the 
geometric onset of the new phase. 

These two descriptions have been shown to be completely equivalent, leading to the same 
critical temperature and the same critical exponents 

Since the percolation approach to the critical behaviour of spin systems is fairly new, we 
briefly sketch its basic features. For simplicity, consider a two-dimensional square lattice 
of linear size R\ on the i?^ sites, we distribute an equal number of up and down spins 
according to the Ising Hamiltonian Hj{T). With decreasing temperature T, adjacent 
sites of parallel spins will begin to form growing clusters, and for T < T^, either the up 
or the down spins will percolate, i.e., form an infinite cluster in the thermodynamic limit 
i? — >■ cx). We thus can define the percolation strength 



which measures the probability that an arbitrary site belongs to the infinite cluster. Since 
P(T) = for all T > and non-zero for all T < T^, it constitutes an order parameter 
for percolation. Although the critical temperature for this transition coincides (in 
two dimensions) with the Curie temperature, the critical exponents in general do not. 
To achieve complete equivalence, including identical critical exponents, the definition of 
cluster has to be modified such that percolation clusters coincide with the correlation 
clusters obtained in the Ising model [0]. This is achieved by assigning to pairs of adjacent 
aligned spins in a geometric cluster an additional bond correlation, randomly distributed 
with the density 



where 2J just corresponds to the energy required for flipping an aligned into a non- 
aligned spin. The modified percolation clusters then consist of aligned spins which are 
bond connected. Only for T = are all aligned spins bonded; for T > 0, some aligned 
spins in a purely geometric cluster are not bonded and hence do not belong to the modified 
cluster. This effectively reduces the size of a given cluster of aligned spins or even cuts 
it into several modified clusters. With such a cluster definition, cluster percolation and 
spontaneous magnetisation become fully equivalent as critical phenomena. 

We now want to extend the percolation definition of critical behaviour to deconfinement 
in SU{2) gauge theory, where we consider clusters of positive (or negative) Polyakov loops 
Li instead of up or down spin clusters; here i denotes the space position of the Polyakov 
loop. First studies have been carried out in two space dimensions 0, using 




T <Tr 



(1) 



rife = 1 — e 



2J/kT 



(2) 



rife = 1 — e 



(3) 
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as bond weight between the (same-sign) Polyakov loops Li and Lj on adjacent spatial 



lattice sites i and j. The effective coupling k is known for small temporal lattices |I0| ; 
we therefore employ for our first study = 2. Using finite size scaling to extrapolate to 
the thermodynamic limit of large systems, our study provides two main results: 

• the percolation coupling ("temperature") coincides with the deconfinement value 
obtained using the expectation value of the Polyakov loop as order parameter; and 

• the critical exponents for percolation strength and cluster size coincide with the 
Ising values predicted by the Svetitsky-Yaffe conjecture for SU{2) thermodynamics. 

For the specific system studied [SU{2) gauge theory, two space dimensions, A^,- = 2], 
we thus find that deconfinement can be specified in terms of Polyakov loop percolation 
with as good or better precision than through the spontaneous symmetry breaking 
measured by Polyakov loop expectation values 0,3. The next steps in this program are 
quite clear - one has to extend Nr to arbitrary values and consider also systems in three 
space dimensions; then the entire study has to be repeated for SU (3) gauge theory. 

One particular advantage of specifying deconfinement in terms of percolation is that this 
appears readily generalizable to full QCD In the presence of dynamical quarks, the 
global Zn symmetry of the pure gauge theory Lagrangian is explicitly broken, so that in 
full QCD the expectation value of the Polyakov loop no longer serves as an order parameter 
for a transition - it is always non- vanishing, even though it continues to vary rapidly in a 
narrow temperature range. In contrast, the per- 
colation strength P{T) remains an order param- 
eter (for the percolation transition) even for fi- b 
nite quark mass niq < oo. This situation corre- ^ 
sponds to the Ising model with non-vanishing ex- 
ternal field if; here there is also no more spon- 
taneous symmetry breaking, with the magnetisa- 
tion m{T, H ^0) 7^ for all T. Nevertheless, the 
percolation transition persists for all values of the 
temperature, tracing a curve Tp[H) of critical tem- 
perature (the so-called 'Kertesz line' [8,12]) in the 
T — H plane. Its counterpart in QCD defines the 
deconfinement diagram shown in Fig. 1, separating 
the confined and the deconfined phase by a line of 
genuine (percolation) critical behaviour, not just 
by a rapid cross-over. It is obviously of great in- 
terest to see if this transition occurs at those tem- 
perature values where the expectation value of the 
Polyakov loop shows its rapid variation. Moreover, 
the relation between the percolation transition and 

the chiral transition at = poses an interesting Figure 1. The Kertesz line of QCD 
open problem. 
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1.2. Colour Screening and Quarkonium Suppression 

In media so dense that several hadrons overlap, the 'other' colour charges screen the 
binding forces between the quarks of any given hadron. As a result, the constituents of the 
medium become deconfined, since screening suppresses the long-range confining potential 

ar ar \ , (4) 

V /ir y 

in a form similar to the Debye screening of the Coulomb potential. The screening mass 
is the inverse of the screening radius; it increases with temperature, so that the range 
of the potential decreases, and eventually deconfinement sets in. Finite temperature 
lattice QCD determines the critical temperature for deconfinement; it also indicates that 
/ic = /u(Tc) ~ 0.4 — 0.6 GeV. With 0.3 - 0.5 fm, the critical colour screening radius is 
therefore considerably smaller than the separation between quarks inside a hadron (about 
1.5 - 2 fm). 

Colour screening provides a probe of the quark-gluon plasma [|l^]. Sufficiently small 
bound states will not melt directly at the deconfinement point; they require a shorter 
screening radius and hence a hotter quark-gluon plasma to be dissociated [0,^. Hence 



quarkonia can be used to to measure the temperature of deconfining media. As exam- 
ple, consider charmonium dissociation. The screening masses for the dissociation of the 
x(3.5) and the J/ip{?>.l) are about 0.35 GeV and 0.70 GeV, respectively ||14[. The x thus 



melts essentially at deconfinement, while the J/ip requires a somewhat higher tempera- 
ture. Similar arguments show that the T melts at a very much higher temperature. For 
quarkonium states in a uniform medium we thus obtain a hierarchy of suppression. 
Summarizing the first section of this survey, we note 

• percolation provides a natural specification and a prospective general order param- 
eter for the onset of deconfinement; 

• the dissociation of quarkonia by colour screening signals deconfinement and serves 
as thermometer of the deconfined medium. 

In the following section, we turn to percolation in nuclear collisions and to the resulting 
suppression of x and direct J/ip production. 



2. PARTON PERCOLATION[|T 

In a high energy nuclear collision viewed in the overall center of mass, the two Lorentz- 
contracted nuclei quickly pass through each other. After about 1 fm at the SPS and only 
about 0.1 fm at RHIC, they have separated and left behind the partonic medium which 
is our potential QGP candidate. Inside this medium, primary nucleon-nucleon collisions 
have left probes, such as heavy quark pairs or Drell-Yan dileptons, whose fate can provide 
information about the state of the partonic medium. 

The average number of produced partons will depend on the number of nucleon-nucleon 
collisions and/or the number of participating ('wounded') nucleons. Up to SPS energies. 



the number of wounded nucleons appears to be the main determining factor [0]. This 



is understandable in terms of the dilated and hence large soft parton formation time 
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as seen in the rest frame of either nucleus. Interference and cancellation effects of the 
Landau- Pomeranchuk type thus prevent soft parton emission at each nucleon-nucleon 
collision. At higher energies, large additional contributions can come from hard partons 
(minijets or jets) with short formation time, and these will be proportional to the number 
of collisions. Even at SPS energy, nucleon stopping and secondary multiplicities at mid- 
rapidity seem to increase somewhat with increasing mass number A, and this could be the 
onset of collision dependent effects. The s-dependence of the mid-rapidity multiplicity is 
presumably a combination of the s-dependence of the parton distribution function and the 
onset of significant hard hadron production. We shall try to include these dependences 
here in rather phenomenological terms, without specifying their origin. 

The number of gluons of transverse momentum fc^ emitted by a wounded nucleon can 
be estimated using the nucleonic gluon distribution function g{x) from deep inelastic 
scattering studies, 

where we have assumed factorization of the parton distribution in x and kx, and taken 
f{k^) as normalized to unity. At ?/ = 0, x ~ kTl\fsi where ^/s is the incident collision 



energy. Using the MRS-H form of the gluon distribution function and integrating the 
resulting Eq. over fc-r then gives us {dN / dy)y=Q ~ 2 for SPS and 4 for RHIC energy. 

Each parton of transverse momentum kx has an effective transverse size r ~ fc^^, 
and nucleus-nucleus collisions provide many such partons overlapping in the transverse 
plane. We want to study their percolation behaviour. To simplify matters, we shall con- 
sider 'average' partons (rather than averaging results with f{k^)). From the transverse 
momentum dependence of Drell-Yan dilepton production through quark-antiquark anni- 
hilation, or from that of J/ ip production through gluon fusion, we know that the effective 
intrinsic transverse momentum of partons is {kx) — 0.75 GeV/c, leading to an average 
transverse parton radius r ~ 0.27 fm. We now consider the percolation pattern for discs 
of that size. 

As a prelude, we assume N discs of radius r to be randomly distributed on a flat surface 
of radius R»r and study the average cluster density rid^n) and the average cluster size 
Sc\{n) as a function of the overall density n = N/ttR^. The result is shown in Fig. 0, 
where we have normalized all quantities to the inverse disc area l/vrr^ in order to make 
them dimensionless. We note that the derivative of the cluster density peaks sharply at a 
certain density; in the 'thermodynamic' limit R — > cxd, it diverges at the percolation point 
Up = 1.175 We see in particular that if we want to reach a certain cluster density 
Uci, attained at a certain overall n, then at this point the average cluster Sci{n) also has 
a certain finite size. In other words, requiring a specific density for the onset of a new 
phase implies that this onset occurs for a specific finite fractional size of the system. The 
percolation point, for example, is reached when about 50% of the surface is covered by 
discs. 

Turning now to nuclear collisions, we have to distribute the partonic discs not on 
a flat surface, but instead in the transverse plane according to the nucleon distribution 



determined by the profiles of the colliding nuclei ||2T|. Moreover, we have to allow collisions 



at different impact parameters, and we want to consider different A — B collisions. 
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Figure 2. Average cluster density nci(n) (left) and average fractional cluster size 
S'ci (n) /S'totai (right) as function of the overall density n of discs, for r/R = 1/20; in 
(b), the derivative of 5'ci(n)/5'totai with respect to n is also shown (dotted line). The 
percolation point in the limit r/R ^ is indicated by Up. 




Figure 3. Cluster density vs. parton den- 
sity for different centralities and different 
A — B configurations at SPS energy. 



Figure 4. Fractional cluster size vs. parton 
density, together with its derivative, for 
Pb—Pb collisions at SPS energy. 
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Figure 5. Cluster density vs. parton den- 
sity for different centralities and different 
A — B configurations at RHIC energy. 
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Figure 6. Fractional cluster size vs. parton 
density, together with its derivative, for 
Cu — Cu collisions at RHIC energy. 



In Fig. ^, the cluster density rid is shown as function of the overall density nparton 
of 'average' partons of radius r = 0.27 fm, for different centralities of various A — B 
combinations, all at SPS energy, where we take 2 partons per wounded nucleon. For 
each A — B, the point at highest riparton corresponds to central (impact parameter 6 = 0) 
collisions, the one at lowest nparton to the most peripheral collisions possible. We note 
that by varying A — B and centrality we produce an essentially universal cluster density 
^ci ('^parton)- Now we assume that deconfinement occurs at the percolation point, which 
can be determined by studying the average fractional cluster size Sd/Stotui- In Fig- 
we see the result for Pb-Pb collisions. At the density for which the derivative peaks, 
Nparton — 4.2 fm~^, the cluster density reaches its critical value rid — 6 fm~^. The same 
critical value is found for Sn — Sn and U — U collisions at SPS energy, even though this 
value is attained at quite different centrality in the different configurations (6 ~ 4.5 fm 
for Sn — Sn, 8.5 fm for Pb-Pb and 10 fm for U — U collisions). 

In Fig's. ^ and ^ we show the result of this procedure for RHIC energy, where we have 4 
gluons per wounded nucleon. Again the critical cluster density at deconfinement is found 
to be nci ~ 6 fm~^. Here we have to consider Cu — Cu collisions in order to study the 
onset of deconfinement, since for Pb — Pb all possible centralities are above the threshold. 

We thus conclude that the critical cluster density is, as expected, a universal quantity, 
independent of the choice of A—B, of centrality, and of incident energy. The fractional size 
of the deconfined cluster at threshold does, however, depend on the collision configuration: 
for Pb-Pb at SPS, we have S'd/^totai ^ 0.45, while Cu-Cu at RHIC gives 0.70. In other 
words, at the higher RHIC energy, deconfinement sets in for a larger bubble than at the 
lower SPS energy. 

We further note that at SPS energy, S — S and S — U are below the deconfinement 
threshold even for the most central collisions, while at RHIC even the most peripheral 
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Figure 7. The J/ip survival probability as function of the parton density for Pb — Pb 
collisions at SPS energy; the dotted line corresponds to a medium of uniform parton 
density, the solid line to collisions with parton densities determined by the profiles of the 
colliding nuclei. 
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Pb~Pb collisions are above the threshold. To study the threshold, we thus require heavy 
nuclei (such as Pb) at the SPS, lighter nuclei (such as Cu) or lower incident energy at 
RHIC. 

Having specified the onset of deconfinement through percolation, we now turn to de- 
confinement signatures. As noted above, the x is dissociated essentially at deconfinement, 
so when the cluster density reaches nj;['* ^ 6 fm~^, all x's which are inside the percolating 
cluster disappear. The amount of x suppression thus depends on the fractional cluster 
size at deconfinement, which is about 45 % for Pb — Pb at the SPS, about 70 % for 
Cu - Cu at RHIC and effectively 100 % for Pb - Pb at RHIC. At deconfinement, the 
corresponding fractions of J/ip production through intermediate x states should thus be 
suppressed. 

As already noted, the dissociation of directly produced J/ip^s requires a a larger screen- 
ing mass and thus higher density than available at the point of deconfinement. On a 
microscopic level, J/ip break-up becomes possible only for harder than the average glu- 
ons at deconfinement. In principle, the relevant quantities can be determined by lattice 
QCD studies; however, this requires computer capabilities which are being reached only 
now. For the time, we therefore choose rj/^ = 0.22 fm and ridlJ/ip) = 7.8 fm~^; this 
corresponds to a ratio of transverse dissociation energy densities eri^J / ^j) / ^rix) — 1-6, 



in accord in studies combining potential theory and lattice results |ll^J15| . With these 
parameters, one obtains 20% direct J/ip suppression in Pb — Pb collisions at SPS energy, 
once the required cluster density is reached. The corresponding fractions for RHIC are 
30% for Cu — Cu and 35% for Pb — Pb collisions. 

Schematically, one thus obtains a suppression pattern as illustrated in Fig. |^, where it 
is assumed that 40% of the observed J/ip's come from x decay and 60% are produced 
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Figure 8. The J/ip survival probability as function of the parton density for Pb—Pb (left) 
and for Cu—Cu (right) collisions at RHIC energy; the dotted line corresponds to a medium 
of uniform parton density, the solid line to collisions with parton densities determined by 
the profiles of the colliding nuclei. 



directly. For clarity purposes, we ignore here the small fraction (about 8%) coming from 
if)' decay. The dotted curve corresponds to the suppression which would take place in 
a uniform medium of precisely specified parton density: at the deconfinement point, all 
x's are dissociated, so that the corresponding fraction of decay J/ip's is gone; when the 
density for direct J/ip melting is reached, these disappear as well, leading to complete 
J /if) suppression. In actual nuclear collisions, the medium is not uniform, with denser 
'inner' and less dense 'outer' regions in the transverse plane. Suppression now occurs 
only in the fraction S^/Stotai for the x P^-rt and Sj/^/ Stotai for the direct J/ip part, where 
and Sj/^ denote the clusters in which the respective dissociation density is reached. 
Since these fractions depend on energy and nuclear geometry, the suppression curves 
are different for different experimental configurations. In Fig. |^, the result is shown for 
Pb—Pb collisions at the SPS, from impact parameter 6 = 11 fm to 6 = fm. The difference 
between this curve and the one for a uniform medium thus reflects the surviving J/ijj^s 
produced in the less dense outer regions. 

In Fig. H, similar calculations are shown for Pb—Pb and Cu~Cu at RHIC energy. Here 
we note in particular that Pb — Pb collisions are for all meaningful centralities (6 < 11 
fm) above both the x ^^id the direct J/ip threshold, so that we get a smooth anomalous 
suppression increasing from about 60% at 6 = 11 fm to about 90% at 6 = 0. Combining 
this suppression with the 'normal' pre-resonance absorption in nuclear matter, we thus 
predict for central Pb — Pb collisions at RHIC a, J/ip production rate of less than 5% the 
corresponding unsuppressed rate (excluding possible B decay contributions). It should 
also be noted that this result is based on twice the number of gluons per wounded nucleon 
at RHIC, compared to the SPS value. A larger increase, based on a possible larger hadron 
multiplicity at RHIC, would lead to more J/V' suppression. Similarly, we assume an 



10 



—3 

CO 




20 40 60 80 100 120 



Et [GeV] 

Figure 9. The J/ip survival probability as function of the measured transverse energy Ej^ 
for Pb~Pb collisions at SPS energy; the dashed line corresponds to collisions with fixed 
Ex — b correlation, the solid line includes the experimental Et — b smearing. The circles 
show the 1997, the squares the 1998 data. 



average number {dNg/dy)y=Q of gluons per wounded nucleon. It is conceivable that very 
central collisions reach into the tail part of the multiplicity distribution, with a larger 
number of hadrons and hence also gluons. This would lead to a larger suppression for 
very central collisions. At the SPS, the NA50 collaboration can check if the basis for 
this exists. By combining measurements of hadron multiplicity, transverse energy E^ and 
forward energy -Ezdc? it is possible to study the number of hadrons per wounded nucleon 
as function of E^ and check if there is an increase at highest E^ values. 

Finally we address briefly the comparison of our results on parton percolation to the 
£^r-dependence of the actual data, which contain an additional smearing due to the fact 
that a given Et bin corresponds to a range of impact parameters and hence parton den- 
sities. Including this effect in the standard way |2^, we obtain the result shown in Fig. 
1^. Included in this figure are the 1997 |2^| and the 1998 data with minimum bias 
determined Drell-Yan reference. The high Et points of the 1997 data show an enhance- 
ment due to rescattering effects, which are removed in the 1998 data using a thinner 
target. While reproducing the overall behaviour, our results clearly show deviations from 
the data in detail. As mentioned, the levelling-off of our curve at high Et would have to 
be modified if an increase in the multiplicity per wounded nucleon at high Et should be 
observed. Another possible modification could enter through a density- dependent char- 
monium dissociation. We have here assumed that all x's melt once the critical density 
is reached. Allowing a partial survival chance, which decreases with increasing density, 
would lead to a steeper drop of the suppression with nparton as well as with Et- One 
possible source for such an effect would be the finite life-time of the deconfining medium 
in actual nuclear collision, as seen e.g. in the parton cascade model to be discussed in the 
next section. 
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Summarizing this section, we note that parton percolation provides a consistent frame- 
work to study the onset of deconfinement as well as the onset of charmonium suppression 
as a signature of the transition. It is clear that the compatibility of the parton basis with 
other, soft hadron production processes has to be checked |]18[ . 



3. COLOUR SCREENING IN THE PARTON CASCADE MODEL 

Probing deconfinement in nuclear collisions means studying the early partonic medium 
left behind when the Lorentz-contracted nuclei have passed through each other - a medium 
consisting of partons, of strings, of flux tubes, or of whatever basis is chosen. Any model 
of this type has to 

• define deconfinement, 

• specify the onset of a particular signature such as charmonium suppression, and 

• check that soft hadron production is consistent with the picture. 

In the past section, we had chosen a purely geometric framework, with deconfinement 
defined as percolation of partons of transverse radius r ~ l/^r, and with a given (energy- 
dependent) number of partons emitted from each wounded nucleon. In this section, we 
want to consider a more detailed microscopic model, based on the interaction of partons 



and the resulting evolution of the partonic medium . For this, we shall use the parton 
cascade model but evidently one could also use other approaches, such as HIJING 
27| or the dual parton model pSf . All such models provide a parton or string distribution 



in the transverse plane and thus in principle constitute an alternative to the nuclear profile 
pattern we had used in the previous section. The time evolution of such a distribution 
from the parton cascade model is shown in Fig. and compared to the nuclear profile 



pattern in Fig. |TT]. It is seen there that the parton cascade model in fact provides a 
transverse parton distribution very similar to that based on the distribution of wounded 
nucleons. 

The parton cascade model describes the partonic medium before the onset of confine- 
ment. To determine when confinement sets in for such a medium, we calculate the screen- 
ing mass fi and recall that lattice QCD finds fic — 0.4 — 0.6 GeV at the deconfinement 
point. The screening mass is given by |^ 

/^^ = / d'kVfik), (6) 

where f{k) specifies the momentum distribution of the partons. For a Bose- Einstein 
distribution, with f{k) = (1 - e-''/^)-\ Eq. (|) gives the familiar /x^ = AnasT'^ obtained 
in perturbative QCD. 

Instead, we now insert for f{k) the time-dependent distribution function fpcM(yk,t) 
from the parton cascade model, in order to obtain the non-equilibrium time evolution of 
the screening mass. First preliminary results for central Pb — Pb collisions at RHIC are 
given in Fig. |l2|, indicating the variation of fi as function of the proper medium life-time 
after full nuclear overlap and as function of the transverse radius of the medium. To 
interpret this behaviour, recall that deconfinement as well as x dissociation occurs for 
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Figure 10. Preliminary PCM calculations for the time evolution of the parton population 
in the transverse plane for central Pb—Pb coUisions at the SPS. 
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Figure 11. Preliminary PCM calculations for the radial density of the parton population in 
the transverse plane of central SPS Pb—Pb coUisions at different coUision times, compared 
to that obtained from the wounded nucleon distribution. 
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Figure 12. Preliminary PCM calculations for the time evolution of the colour screening 
mass in different transverse regions of a central Pb — Pb collision at RHIC energy. 

/i > 0.4 - 0.6 GeV, while the dissociation of directly produced J/ip^s requires 0.7 GeV. 
Hence for essentially all production regions, our preliminary results for the screening mass 
exceed the necessary dissociation values for both x and J/ip for some time. In accord 
with the results from parton percolaton, this would seem to indicate almost complete 
J/ijj suppression for central Pb—Pb collisions at RHIC. To really conclude this, however, 
one must first specify how long n has to exceed its critical value in order to 'melt' the 
charmonium state in question. In any case, for really quantitative precision more extensive 
calculations are necessary. 

Let us finally note one particular advantage of microscopic evolution picture such as the 
parton cascade model. Once we have defined the conditions for quarkonium suppression, 
be it through percolation or through colour screening, the parton cascade model can be 
used to check if the conditions needed for the observed suppression lead to consistent 
soft hadron production at the later times. Such models thus provide a bridge between 
the early deconfined and the later hadronic stages. The situation is summarized in Fig. 
0, which also serves as our conclusion. It shows on one hand the general evolution of 
the collision, on the other the effect of this evolution on charmonium production as hard 
probe of the early medium. 
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